Interactions between brown-dwarf binaries and Sun-like 
stars 
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Abstract Several mechanisms have been proposed for 
the formation of brown dwarfs, but there is as yet no 
consensus as to which - if any - are operative in na- 
ture. Any theory of brown dwarf formation must ex- 
plain the observed statistics of brown dwarfs. These 
statistics are limited by selection effects, but they are 
becoming increasingly discriminating. In particular, it 
appears (a) that brown dwarfs that are secondaries to 
Sun-like stars tend to be on wide orbits, a > 100 AU 
(the Brown Dwarf Desert), and (b) that these brown 
dwarfs have a significantly higher chance of being in a 
close (a < 10 AU) binary system with another brown 
dwarf than do brown dwarfs in the field. This then 
raises the issue of whether these brown dwarfs have 
formed in situ, i.e. by fragmentation of a circumstel- 
lar disc; or have formed elsewhere and subsequently 
been captured. We present numerical simulations of the 
purely gravitational interaction between a close brown- 
dwarf binary and a Sun-like star. These simulations 
demonstrate that such interactions have a negligible 
chance (< 0.001) of leading to the close brown-dwarf 
binary being captured by the Sun-like star. Making the 
interactions dissipative by invoking the hydrodynamic 
effects of attendant discs might alter this conclusion. 
However, in order to explain the above statistics, this 
dissipation would have to favour the capture of brown- 
dwarf binaries over single brown-dwarfs, and we present 
arguments why this is unlikely. The simplest inference 
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is that most brown-dwarf binaries - and therefore pos- 
sibly also most single brown dwarfs - form by fragmen- 
tation of circumstellar discs around Sun-like protostars, 
with some of them subsequently being ejected into the 
field. 
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1 Introduction 



numerical methods 



The existence of brown dwarfs was first suggested by 
Kumar (1961), on the basis of theoretical calculations 
of the structure of newly-formed stars, and the realisa- 
tion that stars of sufficiently low mass, as they condense 
out of the interstellar medium, become dense enough 
to be supported by electron degeneracy pressure before 
they become hot enough to burn hydrogen. However, 
it was not until 1995 (Rebolo et al. 1995, Nakajima 
et al. 1995) that brown dwarfs were actually observed. 
Since then brown dwarfs have been observed in increas- 
ing numbers, and in increasing detail, so that we now 
have some idea of their statistical properties - albeit 
subject to strong selection effects due to the fact that 
brown dwarfs are intrinsically faint, and become ever 
fainter with age (Luhman et al. 2007). It is therefore 
appropriate to evaluate how the proposed mechanisms 
for brown dwarf formation measure up against the ob- 
servational constraints. 

In discussing the formation of brown dwarfs, it 
is inappropriate to distinguish them from hydrogen- 
burning stars, since there is no reason to suppose that 
the processes that determine the mass of a star con- 
densing out of interstellar matter change radically at 
the hydrogen-burning limit (~ 0.08 M Q ). Therefore 
brown dwarfs are likely to form a continuum with low- 
mass hydrogen-burning stars. This does not mean that 
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brown dwarfs form in exactly the same way as all other 
stars, but simply that any systematic changes in go- 
ing from massive star formation to brown dwarf forma- 
tion are part of a continuous trend which also takes in 
the formation of intermediate- and low-mass hydrogen- 
burning star^|. Indeed, we propose that, as one goes to 
lower masses, an increasing fraction of stars is formed 
by fragmentation of protostellar accretion discs. 

1.1 Formation mechanisms 

The three main mechanisms for brown dwarf forma- 
tion currently being advocated are the following. Other 
mechanisms have been proposed (Whitworth et al. 
2007), for example the photo-erosion of a pre-existing 
core (Whitworth & Zinnecker, 2005), but are unlikely 
to be major contributors. 

1.1.1 Turbulent generation of very low-mass prestellar 
cores. 

In this mechanism, brown dwarfs can form in isolation 
from exceptionally small and dense prestellar cores. It 
is argued that interstellar turbulence occasionally de- 
livers compressive flows of sufficient strength to form 
such cores (e.g. Padoan & Nordlund 1999; Hcnnebelle 
& Chabrier 2008, 2009). We note that a prestellar core 
of mass M, with isothermal sound speed a, only col- 
lapses if the external pressure exceeds 



x MIN — Q 3 yp 7 W 

For a core at the opacity limit, M ~ 0.003 M©, with 
a ~ 0.2 km s^ 1 (corresponding to molecular gas at T ~ 
10 K), this yiekkS 

^ = 0.13 x 10 10 (H 2 cm- 3 ) (kms- 1 ) 2 ; (2) 

for a core at the brown dwarf limit, it gives 

= 0.20 x 10 7 (H 2 cm- 3 ) (kms" 1 ) 2 . (3) 

m H 2 



1 We note that this is a much less extreme view than that espoused 
by Thies & Kroupa (2007, 2008), who argue that there is an 
abrupt and discontinuous change in the IMF at or just above the 
brown-dwarf limit. 

2 Note that, instead of dividing out the Boltzmann constant to 
give the pressure in terms of a number-density times a tempera- 
ture (as is routinely done), we have divided out the mean mass 
per hydrogen molecule - which for molecular gas with a hydrogen 
mass fraction X = 0.70 is m H2 = 2m p /X ~ 4.8 X 10~ 24 g - so as 
to give the pressure in terms of the number-density of H2 times 
a squared velocity. This highlights the very large H2 number 
densities and velocities required to deliver a ram pressure of this 
order. 



There is presently no observational evidence for the 
dense, supersonic, coherently converging flows required 
to deliver these high ram pressures. Moreover, this 
mechanism has not been simulated, and therefore there 
are no predictions for the statistics of discs, outflows, 
and multiplicity to enable a comparison with observa- 
tion. 

1.1.2 Collapse and fragmentation of larger prestellar 
cores. 

In this mechanism, brown dwarf stars form collec- 
tively along with hydrogen-burning stars, during the 
collapse and fragmentation of more massive prestellar 
cores (i.e. cores which spawn more than one protostar) . 
Brown dwarfs are then delivered to the field by dy- 
namical ejection (Reipurth & Clarke, 2001; Bate, Bon- 
nell & Bromm, 2003; Goodwin, Whitworth & Ward- 
Thompson, 2004a,b; Bate 2009; Offher et al. 2009, 
2010). Only the simulations by Offner et al. (2009) take 
proper account of radiative feedback from newly-formed 
protostars, and, although they form brown dwarfs, (a) 
these stars are still accreting at the end of the simu- 
lation, so their final masses are uncertain, and (b) the 
brown dwarfs formed are too few (i.e. 2) to allow sta- 
tistical (as distinct from phenomenological) inferences 
to be made. It is presently unclear whether this mech- 
anism can deliver the relative number of brown dwarfs 
found in nature (i.e. the stellar IMF), and it offers no 
explanation for the rather particular binary statistics 
of brown dwarfs. 

1.1.3 Fragmentation of protostellar accretion discs. 

In this mechanism, a primary star forms attended by 
an accretion disc. Material infalls onto this disc and 
spirals into the primary protostar at a rate determined 
by how quickly it can redistribute angular momentum. 
However, if the disc becomes sufficiently massive, ex- 
tended and cool, it fragments to produce secondary 
stars. These secondary stars are usually of low mass, i.e. 
very low- mass H-burning stars (VLMSs), brown-dwarf 
stars and planetary- mass stars (planemos) . This mech- 
anism has the merit that it has been simulated numer- 
ically in detail, and reproduces the critical statistical 
properties of brown dwarfs (Whitworth & Stamatellos, 
2006; Stamatellos, Hubber & Whitworth, 2007; Sta- 
matellos & Whitworth, 2009). The VLMSs tend to end 
up closer to the primary (due to a combination of scat- 
tering and secular migration) , most of the planemos and 
brown dwarfs are scattered out of the system. Most of 
the secondaries are attended by small accretion discs - 
large enough to fuel signatures of ongoing accretion, but 
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too small to increase the masses of the secondaries sig- 
nificantly. Many of the secondaries are in close binary 
systems. Many of these close binary systems remain 
on wide orbits round the primary but some are ejected 
and survive the ejection process. Offncr et al. (2010) 
have shown that when radiative feedback is included 
in their simulations, disc fragmentation is suppressed, 
due to the large accretion luminosity from the primary 
protostar. However, Stamatellos, Hubber & Whitworth 
(2011) show that when the episodic nature of accretion 
is taken into account, there are low-luminosity periods 
when disc fragmentation is likely to occur. 

1.2 Observational constraints on the properties of 
brown dwarfs 

1.2.1 The collective properties of brown dwarfs 

There is a growing body of evidence suggesting that 
brown dwarfs form in similar ways to low-mass hy- 
drogen - burning stars. The Initial Mass Functions 
for low-mass hydrogen-burning stars and for brown 
dwarfs appear to merge continuously at the hydrogen- 
burning limit (e.g. Moraux et al. 2007; Lodieu et al. 
2009). Both young brown dwarfs and young low- mass 
hydrogen-burning stars are attended by discs (e.g. Luh- 
man et al. 2008), show signatures of accretion, at rates 
which follow the approximate relationship M cx M 2 
(e.g. Natta et al. 2004; Mohanty, Jayawardhana & 
Basri, 2005; Comcron, Tcsti & Natta 2010), and drive 
jets and outflows (e.g. Whelan et al. 2007). In gen- 
eral, young brown dwarf stars appear to be clustered in 
the same way as hydrogen-burning stars (e.g. Luhman, 
2006). 

1.2.2 The binary statistics of brown dwarfs 

As well as reproducing the collective properties of in- 
dividual brown dwarfs, a theory of brown dwarf for- 
mation must also explain their binary statistics. There 
are three clear trends emerging here, (i) Brown dwarfs 
that are secondaries to Sun-like stars tend to be on 
quite wide orbits (separations S > 100 AU); this is the 
Brown Dwarf Desert (McCarthy & Zuckerman, 2004; 
Sahlmann et al., 2011). (ii) Brown dwarfs that are in 
binary systems with other brown dwarfs (~ 10 — 20% of 
them) tend to be on close orbits (typically S < 20 AU 
and S ~ 5 AU) and to have components of comparable 
mass (q <~ 1) (e.g. Bouy et al., 2003; Burgasser et al., 
2003; Close et al., 2003; Ahmic et al., 2007; Reid et al., 
2008; Gelino & Burgasser, 2010; Todorov, Luhman & 
McLeod, 2010). (iii) A brown dwarf that is a secondary 
in a wide orbit around a Sun-like star is estimated to 
be at least twice as likely to be in a close binary with 



another brown dwarf (and - approximately - five times 
as likely to be in a triple, and thirteen times as likely 
to be in a quadruple, with other brown dwarfs) as a 
brown dwarf in the field (e.g. Burgasser, Kirkpatrick 
& Lowrance, 2005; Faherty et al., 2010; Faherty et al., 
2011). As long as these statistical inferences hold up in 
the face of improved observational data and analysis, 
and larger sample sizes, they constitute very discrimi- 
nating constraints, which theories of the formation of 
brown dwarfs should seek to address. 

The fact that the observations are biased towards 
younger brown dwarfs (because younger brown dwarfs 
are brighter), is not a critical issue in this context, 
since it weights the statistics towards those properties 
of brown dwarfs that reflect how they are formed, as 
distinct from those that are produced by subsequent 
longer-term environmental and evolutionary effects. 

1.3 The origin of brown-dwarf binaries 

If brown dwarfs form from very low-mass prestel- 
lar cores, close brown-dwarf binaries with appropriate 
masses and separations may form by secondary frag- 
mentation when molecular hydrogen dissociates, as ar- 
gued by Whitworth & Stamatellos (2006). However, 
this mechanism has never been demonstrated to work 
in numerical simulations (e.g. Bate, 1998). Moreover, 
if it does work in nature, it only delivers isolated brown- 
dwarf binaries, and there remains a problem with ex- 
plaining the brown-dwarf binaries that are found in 
wide orbits round Sun-like stars. The results we re- 
port here suggest that, at least in the absence of hydro- 
dynamic dissipation, the capture of close brown-dwarf 
binaries into wide orbits round Sun-like stars is very 
unlikely. 

If a brown dwarf forms along with other stars from a 
single core, the subsequent evolution is highly chaotic, 
and a range of outcomes is possible. The result- 
ing stars tend to define themselves (in terms of their 
masses) during the highly dissipative phase when the 
gas switches from isothcrmality to adiabaticity at den- 
sities p ~ 10~ 13 gcm~ 3 . These high densities may not 
arise exclusively near the centre of mass of the core, 
but also in filaments channeling material into the cen- 
tre (e.g. Walch et al. 2010; Girichidis et al. 2010), and 
it is unclear where conditions will be most favourable 
for forming brown dwarfs. However, whether they form 
near the centre (where the densities are higher) or in 
a filament away from the centre (where the temper- 
atures are lower), they quickly end up in the centre. 
This paper demonstrates that, if the violent interac- 
tions with other stars that then occur are purely gravi- 
tational, they are unlikely to lead to the formation and 
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survival of brown-dwarf binaries, and even less likely 
to deliver such binaries into wide orbits around more 
massive stars. 

If brown dwarfs form by disc fragmentation, their 
statistical properties appear to meet all the critical 
observational constraints (Stamatellos & Whitworth 
2009). 

• Discs provide an ideal place to form objects close to 
the opacity limit. The requirement that the cooling 
timescale in a fragment be less than the dynamical 
timescale (Gammie, 2001), which is, in effect, the 
opacity limit, is more easily met in a disc, since ma- 
terial is parked in the disc and can cool quasistatically 
until it becomes Toomre unstable (Toomre 1964). By 
contrast, in a violently compressive turbulent flow, 
the dynamical timescale, and hence the time avail- 
able for cooling, are by definition shorter, and so an 
adiabatic bounce is more likely. 

• Disc fragmentation naturally produces protostars 
with low masses and a mass-function dM /dM oc 
M-°- 6 . 

• The distribution of angular momentum in a disc en- 
sures that the reservoir of material that a particular 
protostar can accrete is limited, so there is much less 
need to cut off accretion by dynamical ejection. Ejec- 
tions do occur, but there is no need for them to occur 
rapidly. 

• Fragmentation only occurs in the outer reaches of a 
disc (R > 100 AU). The VLMSs that form by frag- 
mentation tend to be scattered or migrate into closer 
orbits, but the brown dwarfs and planemos tend to 
be scattered outwards and/or ejected. 

• Protostars formed by disc fragmentation are at- 
tended by their own small accretion discs, with 
masses and radii comparable with those inferred from 
observation. 

• Protostars condensing out of discs have a natural ten- 
dency to pair up and form close binaries, and these 
close binaries tend to be on wide orbits around the 
primary (Sun-like) star. 

• Many of the brown dwarfs and planemos are dynam- 
ically ejected. However, since these ejections involve 
rather mild interactions with other low-mass frag- 
ments, accretion discs and close low-mass compan- 
ions can survive the ejection process. 

It has been argued (a) that when radiative feedback 
from the central primary protostar is included, disc 
fragmentation is suppressed (Offher et al. 2009, 2010), 
and (b) that fragmenting discs are not seen (Maury et 
al. 2010). It is therefore appropriate to comment briefly 
on these issues, (a) Stamatellos, Hubber & Whitworth 
(2011) show that if accretion onto the central primary 



protostar is episodic - and there is considerable evi- 
dence to suggest that this is the case - then there are 
periods during which the accretion rate, and hence the 
luminosity, are low, and these periods are sufficiently 
long to allow the disc to become very cool and frag- 
ment, (b) The timescale for dynamical fragmentation is 
very short for brown dwarfs, i DYN ~ 6 kyr (M/0.01M Q ). 
Therefore the probability of observing a fragmenting 
disc with current technology is very low - although it 
will improve with Herschel and ALMA surveys (Sta- 
matellos et al. 2011). 

In the meantime, it is appropriate to explore the pos- 
sibility of forming brown-dwarf binaries other than by 
disc fragmentation, and then capturing them into wide 
orbits around Sun-like stars. This paper reports nu- 
merical experiments which suggest that this is unlikely. 
In Section [2] we define the initial conditions used for the 
numerical experiments. In Section[3]we describe the nu- 
merical method. In Section |4] we define the categories 
used in interpreting the results. Section [5] presents the 
results, and Section [5] presents discussion of the results 
and summarises our conclusions. 



2 Initial conditions 

We follow the interaction of three point-mass stars, la- 
belled 1, 2, 3. The first star is a Sun-like star with mass 
Afj = 1 M Q , and the other two are brown dwarfs with 
masses M 2 = M 3 = 0.05 M©. 

At the start of a simulation, the two brown dwarfs 
(Stars 2 and 3) are in a close binary system, on a cir- 
cular orbit with separation S 23 = 5 AU, hence period 
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2n 
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2GAL 
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35 years , 



orbital velocity 



v„ = I 



GM 2 
and net energy 



1/2 
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GM 2 M 3 



2.1kms~ 



■ 4.6 x 10 42 erg . 



(4) 



(5) 



(6) 



These are representative parameters for a browm-dwarf 
binary in the field. The centre of mass of this binary 
is placed at the origin of coordinates. The orientation 
and phase of the binary orbit are random. 

At the start of a simulation, the Sun- like star (Star 1) 
has position (x,y,z) = (b, — 3000, 0) AU with b chosen 
randomly from 



Pb db 



2bdb 



(3000 AU)' 



0AU < b < 3000AU;(7) 
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and velocity {v x ,v y ,v z ) — (0, |v|,0)kms 1 with v cho- 
sen randomly from 



p v dv 



exp I — — — )■ (//' . 



(2tt) 1 /2 (T 



2<t 2 



(8) 



with o~ = 1, 2, 3, 4kms l . Hence the additional system 
energy introduced by the Sun-like star is 



E, 



GM X (M 2 +M a 



M 1 (M 2 +M 3 )w 2 



((3000 AU) 2 + 6 2 ) 1/2 2 (M, + M 2 + M 3 ) 

2 X -1/2 



0.3 x 10 42 erg I 1 



3000 AU 



°- 9Xl ° 42erg (kn7^) : 



(9) 



The velocity dispersions are representative of those in 
young star clusters, and we are assuming that the ve- 
locities are random and isotropic, with a Gaussian dis- 
tribution. The distribution of impact parameters takes 
account of the fact that the integrated cross-section at 
impact parameter b is irb 2 ; the range of impact param- 
eters is that for which significant interactions between 
the Sun-like star and the brown dwarf binary can be ex- 
pected. The starting positions are chosen so that grav- 
itational focussing is captured. We note that, although 
the approach of the Sun-like star is in the z = plane, 
these are not coplanar interactions, since the orienta- 
tion of the orbit of the brown dwarf binary is random. 

In our numerical experiments, we only vary five pa- 
rameters: two angles, (#,</>), giving the orientation of 
the orbit of the brown-dwarf binary; a third angle, ip, 
giving the phase of one of the brown dwarfs in its or- 
bit; and the impact parameter, 6, and velocity, v, of 
the incoming Sun-like star. We do not consider differ- 
ent masses for the stars (Mi, M2, -M3), nor do we con- 
sider close brown-dwarfs binaries with different initial 
orbital parameters (separations, S23, and eccentricities, 
623). The reason for limiting our numerical experiments 
to the first five parameters is that the outcome of an 
interaction depends in a chaotic manner on the orien- 
tation of the brown-dwarf orbit, the phase of the brown 
dwarfs in this orbit at the point of closest approach of 
the incoming Sun-like star, and the velocity and dis- 
tance of the incoming Sun-like star at the point of clos- 
est approach. If one of the last five parameters were 
changed, the areas of this phase-space corresponding to 
different generic outcomes (see Section [5] and Table [lj 
would shift, but the overall topology of the phase-space 
and the relative frequencies of different outcomes would 
change very little. In other words, if the Sun-like star 
considered here almost never captures the close brown- 
dwarf binary considered here, it is extremely unlikely 



that Sun-like stars of different (but comparable) mass 
would capture brown-dwarf binaries with different (but 
comparable) properties. 

We reiterate that the collision parameters are chosen 
to reflect the typical velocity dispersions in star clusters 
(lkms -1 < a < 4kms _1 ) and the range of impact pa- 
rameters (0 < b < 3000 AU) for which a close interac- 
tion is possible, but not guaranteed, at these velocities. 
In effect, this means exploring impact parameters and 
starting distances at which the gravitational potential 
energy of the brown-dwarf binary relative to the Sun- 
like star is less than their relative kinetic energy. Under 
this circumstance, there are many collisions that are 
simply flybys, and leave the brown-dwarf binary intact 
(outcome S in Table [T]). Obviously our main concern is 
with the non-flybys. 

For each value of a we perform a Monte-Carlo en- 
semble of 10 6 realisations, each of which is followed for 
150 kyr. A single realisation is defined by the five ran- 
dom numbers used to generate 0, <fi, tp, v and b. 



3 Numerical method 

Each realisation is evolved using a fourth-order Hermite 
scheme with an adaptive global timestep given by 



At = 7 MIN 



IAv. 



|Ar. 
|Aa. 



1/2' 



(10) 



here Ar. . = r r , Av. . = v v , and Aa. . = a a , 
and the minimisation is taken over all pairs of stars, 
ij = 12, 23, and 31. 

We use 7 = 0.01 , and this ensures that the frac- 
tional error in the total energy, E TOT , linear momen- 
tum, P TOT , and angular momentum, H TOT , satisfy 



log 
log 

°io 



log 



|A£L 



I -^"tot I 

|AP TOT 
IP 

I" 1 - TOT I 

|AH T ,„ 



|H„ 



-5.5 ±1.0. 



-13.3 ± 1.2. 



-8.6 ± 1.5. 



(11) 
(12) 

(13) 



These uncertainties can be reduced further, but this 
would not significantly alter the statistics of the en- 
semble. 



4 Definitions 

At any time, and for each pair, ij, we can compute the 
total mass, M^-, centre-of-mass position, ry, centre-of- 
mass velocity, Vy and mutual energy, according 
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Fig. 1 The left hand figures give the numbers of systems in the different designations (X, NHT, HT2, HT1 and D), as 
a function of time; we omit designation S, since, although it is the most common outcome, it corresponds to cases where 
there is no significant interaction between the Sun-like star and the brown-dwarf binary. These plots show that at the 
end the dominant designation after S is X (one brown dwarf captured by the Sun-like star, one ejected); followed by NHT 
(non- hierarchical triple); followed by HT2 (hierarchical triple with one brown dwarf in a close orbit round the Sun-like star, 
and one in a wide orbit). There are very few HT1 systems (hierarchical triples with the brown dwarfs in a close system 
on a wide orbit round the Sun-like star) or D systems (total disruption). The right hand figures give the rates at which 
systems are exchanged between designations. This shows that the main channels are converting NHT systems into HT2 or 
X, and HT2 into X. In other words, in almost all cases where there is a significant interaction, one of the brown dwarfs 
ends up either ejected directly , or scattered into a wide orbit and then ejected, or scattered into a wide orbit. The first 
row presents results for a v — lkms -1 , the second row for a v = 2kms~ , the third row for a v — 3kms _1 , and the fourth 
row for <7„ = 4kms -1 . 
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Table 1 Abbreviations for different outcomes, and their definitions. 



Abbreviation 


Definition 


HT1 


Hierarchical Triple of Type 1 




(close system contains both brown dwarfs); 


HT2 


Hierarchical Triple of Type 2 




(close system contains the Sun-like star); 


NHT 


Non-Hierarchical Triple; 


X 


eXchange (one brown dwarf captured 




by Sun-like star, other ejected); 


S 


Survival (brown-dwarf binary survives, 




but not captured by Sun- like star); 


D 


Disintegration (all stars become single, 




only possible with extreme v). 



to 



My = Mi + Mj , 

+ M,v, 



MjVj + Mm 



E 



GMiMj MjMjlvi-Vj-1 



2 (Mi + Mj) 



(14) 
(15) 

(16) 

(17) 



We can also compute the energy of the third star, k, 
relative to this pair (treated, rightly or wrongly, as an 
unresolved pair): 



E k = 
If 

E k > 0, 



GMijMk , MijM k \ 



2 (My + M k ) 



(18) 



(19) 



we label k as an escaper (but see below). If k = 1, 
we designate the system S, meaning that the Sun- like 
star has remained single, and the brown-dwarf binary 
has survived intact (albeit it with an adjusted orbit). If 
k =/= I, we designate the system X, meaning that there 
has been an exchange, one of the brown dwarfs has been 
ejected and the other has been captured by the Sun-like 
star. If none of the E k satisfies Inequality (fl"9)). but one 
of the pairs satisfies 



Sun-like star is presumed to have captured the brown- 
dwarf binary into a (relatively) wide orbit; in HT2 sys- 
tems, the close pair is ij = 12 or ij = 13, and therefore 
the Sun- like star has captured one of the brown dwarfs 
into a close orbit, and scattered the other one into a 
wide orbit. If neither Inequality (|19[) nor Inequality 
(|20|) are satisfied, the system is designated NHT for 
non-hierarchical triple. We summarise these categories 
in Table [Q 

These labels and designations are not rigorous, in 
the following sense. There are transient configurations 
in which star k satisfies Inequality (fl~9"|) but is not really 
an escaper, or pair ij satisfies Inequality (|20|) but the 
system in not really hierarchical. For example, one of 
the stars in a close binary system may have such a large 
orbital velocity, and hence such a large velocity relative 
to the third star, that it appears to be an escaper, even 
though the binary as a whole is bound to the third star, 
and vice versa. However, these circumstances are both 
short-lived and increasingly rare, so that once the sys- 
tem has been evolved for 150 kyr (i.e. more than 4000 
times the initial binary period), they have a negligible 
impact on the statistics. 



5 Results 

Fig. [1] presents the results obtained with a = lkms -1 
(first row), 2kms -1 (second row), 3kms _1 (third row) 



Ea < 10MW{E ik ,E jk ) , 



(20) 



we label ij as a close pair (but again see below), and 
designate the system hierarchical. The factor 10 is 
somewhat arbitrary, but the results do not depend 
strongly on the number chosen here; the larger the num- 
ber, the more stable are the hierarchical systems. We 
distinguish two types of hierarchical system: In HT1 
systems, the close pair is ij = 23, and therefore the 



Table 2 Relative frequencies of different outcomes. 



a v /(km/s) 


1 


2 


3 


4 


/s 


0.7069 


0.8243 


0.8726 


0.8994 


fx 


0.2564 


0.1541 


0.1129 


0.0895 


f NHT 


0.0249 


0.0143 


0.0097 


0.0074 


/ HT2 


0.0117 


0.0066 


0.0044 


0.0035 


/hti 


0.0009 


0.0005 


0.0003 


0.0002 


f D 


0.0002 


0.0002 


0.0001 


0.0000 
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and 4kms~ 1 (fourth row). On each row, the left panel 
shows the number of systems in each category apart 
from S (see Table H}, as a function of time, and the 
right panel shows the rates at which systems are trans- 
fered from one category to another (only the dominant 
channels, i.e. NHT^X; NHT— »HT2; HT2^X), again 
as a function of time. 

Table [5] summarises the final outcome in terms of 
the fraction of experiments that produce systems in the 
different categories after 150,000 years. The dominant 
category is S (survival, i.e. the brown-dwarf binary sur- 
vives and is not captured by the Sun-like star) , followed 
by X (exchange, i.e. one brown dwarf is captured by the 
Sun- like star and the other is ejected), followed by NHT 
(non-hierarchical triple), followed by HT2 (hierarchical 
triple of type 2, i.e. one brown dwarf captured into 
a close orbit around the Sun-like star, and the second 
brown dwarf captured into a wide orbit around the Sun- 
like star), followed by HT1 (hierarchical triple of type 
1, i.e. the close brown-dwarf binary survives as such 
but is captured into a wide orbit around the Sun-like 
star), followed by D (total disruption, i.e. three single 
stars). Thus the majority of systems that experience a 
significant interaction undergo exchange and end up as 
X systems, i.e. one of the brown dwarfs is captured by 
the Sun-like star, and the other is ejected. There is an 
extremely low probability of producing HT1 systems, 
i.e. systems like the observed ones with a close brown- 
dwarf binary on a wide orbit around a Sun-like star. 
Furthermore, for each such system produced, there are 
many more HT2 and X systems produced, i.e. hier- 
archical systems with one brown dwarf in a close orbit 
around the Sun-like star and the other on a much wider 
orbit, and exchange systems where one of the brown 
dwarfs has been ejected. 

As the velocity dispersion, a v , is increased, the 
approach velocity of the Sun-like star increases, and 
the influence of gravitational focussing diminishes, so 
there are fewer significant interactions. Consequently 
df s /da v > 0, but for all other outcomes (O = X, NHT, 
HT1, HT2, D) df Q /da v < 0. 

6 Discussion and conclusions 

The numerical experiments described above suggest 
that there is an extremely low probability that a Sun- 
like star captures a pre-existing close brown-dwarf bi- 
nary into a wide orbit, if the interaction between them 
is purely gtavitational. If there is such an interaction, 
the end result is almost always that one brown dwarf is 
cither ejected, or scattered into a wide orbit, whilst the 
other becomes a close companion to the Sun-like star. 



One might argue that, if the Sun-like star and/or 
the brown-dwarf binary were attended by circumstellar 
and/or circumbinary discs, interactions between them 
would be more dissipative and it might then be possi- 
ble for the Sun-like star to capture the brown-dwarf 
binary into a wide orbit. However, disc dissipation 
would require a rather close interaction, and this would 
be likely to leave the brown-dwarf binary on an ec- 
centric orbit with a close periastron to the Sun-like 
star. With the disc now dissipated, the simulations pre- 
sented here indicate that it is then virtually inevitable 
that one of the brown dwarfs will end up either being 
ejected, or being scattered into a wide orbit, whilst the 
other becomes a close companion to the Sun-like star. 
This would be in direct conflict with the observed lack 
of brown dwarfs orbiting close to Sun-like stars (The 
Brown Dwarf Desert). 

Furthermore, if close brown-dwarf binaries retain 
dissipative discs long enough for them to regulate close 
interactions with Sun-like stars, it seems likely that sin- 
gle brown dwarfs do too. However, this will make it 
easy for single brown dwarfs to be captured into close 
orbits around Sun-like stars, and again this will be in 
conflict with The Brown Dwarf Desert. 

Another possibility is that a brown dwarf binary 
might be delivered into a wide orbit around a Sun-like 
star by an exchange reaction, i.e. a Sun-like star which 
already had some other sort of companion might swap 
this companion for a brown-dwarf binary in a four-body 
interaction. We cannot comment quantitatively on this 
possibility on the basis of the results presented here 
(which only involve three bodies), but it seems rather 
contrived and unlikely. First, such exchanges normally 
work in the direction of swapping a companion having 
a low mass for one having a higher mass; therefore this 
route requires some other mechanism to form the ini- 
tial binary system with a Sun-like primary and a very 
low-mass secondary, as a prelude to swapping the very 
low-mass secondary for a brown dwarf binary. Second, 
in an environment where exchanges occur (i) such ex- 
changes will also act to remove a brown dwarf binary 
bound to a Sun-like star, in favour of an even more 
massive companion, and (ii) non-exchange interactions 
will scatter brown dwarf binaries that remain bound to 
Sun-like primaries into a range of orbits, thereby pop- 
ulating The Brown Dwarf Desert. 

We conclude that close brown-dwarf binaries in wide 
orbits around Sun-like stars are likely to have been 
formed there, by disc fragmentation, and that if these 
brown dwarfs formed by disc fragmentation, then other 
brown dwarfs (those not in close binary systems and / or 
no longer in orbit around Sun-like stars) may also have 
formed by disc fragmentation. We note that Simula- 
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tions of disc fragmentation reproduce all the discrimi- 
nating statistical properties of brown dwarfs (Stamatel- 
los & Whitworth 2009), and we advocate that alterna- 
tive theories of brown dwarf formation be simulated 
with sufficient resolution to enable the same compar- 
isons to be made. 
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